Anisotropic critical fields of MgB2 single crystals * 
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The recently discovered superconductivity in MgB2 has created the 
world sensation. In spite of the relatively high superconducting transition 
temperature T c = 39 K, the superconductivity is understood in terms of 
rare two gap superconductor with energy gaps attached to the a and ir- 
band. However, this simple model cannot describe the temperature depen- 
dent anisotropy in H C 2 or the temperature dependence of the anisotropic 
magnetic penetration depth. Here we propose a model with two anisotropic 
energy gaps with different shapes. Indeed the present model describes a 
number of pecularities of MgB 2 which have been revealed only recently 
through single crystal MgB 2 . 

PACS numbers: 74.20.Rp, 74.25.Bt, 74.70.Ad 
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1. Introduction 

The discovery of new superconductivity in MgB2 took the world by 
surprise [1]. Early studies based on polycrystalline samples lead to a two 
gap model [2, 3, 4]. On the other hand the anisotropy in H C 2(t, G) suggested 
an anisotropic s-wave model [5, 6, 7, 8, 9]. Further it is clear that the 
simple two gap model cannot describe the strong temperature dependence 
of j(t) = H$(t)/HZ 2 (t) observed in single crystal MgE>2 [6, 9]. For this 
we need an order parameter A(/c) of oblate shape [9]. On the other hand, 
it is well known that H$(t)/H£ 2 (t) > 1 and H^(t)/H^(t) > 1 for single 
crystal experiments [10, 11, 12], which contradicts the Ginzburg-Landau 
phenomenology. Further, both magnetic penetration depth data and H^(t) 
data suggest a prolate order parameter as in [7, 8]. Indeed, an earlier STM 
study suggested a prolate order parameter as well [13]. 

Is the order parameter prolate or oblate? The answer is that we need 
both. We suggest, that the oblate order parameter is attached to the cr-band 
while the prolate order parameter to the 7r-band. In the following we shall 
describe salient features of single crystal MgE>2 within the present model. 

2. Upper critical field 

We just point out two features in the upper critical field which are 
outside of the Ginzburg-Landau phenomenology. 

(a) the strong temperature dependence of the anisotropy parameter 
7 (t) = H$(t)/H£ 2 (t) [9]. (b) second the deviat ion from the effective mass 
model H c2 (t,@)/H c2 (t,0) + (cos 2 G + 7 sin 2 G) -1 / 2 [6]. 

For example A (A;) ~ l/\/l + az 2 , where z = cosG and a ~ 100 can 
describe the temperature dependence of 7(t), as has been shown in Ref. [9]. 
Also, by fitting the experimental data we obtained v a ~ 2.7 • 10 7 cm/sec 
and v c /v a ~ 0.48. We point out that these values are very consistent with 
the ones for the cr-band [2]. 

By fitting experimental data for H C 2{6) [14] close to T c we can deduce the 
ratio v c /v a ~ 0.73, implying a = 10, for example. This result is bigger than 
that estimated in our previous analysis of H C 2(t) [9]. Next order correction 
in (T — T c )/T c leads to the increase of the ratio v c /v a for the same values 
of the parameter a. In the future we are going to elaborate our study of 
the upper critical field behavior in MgB2 by taking into account impurity 
scattering, what we believe will improve the agreement with experiment. 

3. Lower critical field 

^,From c-axis oriented films and single crystals, the superfluid density 
p s ^ a and the lower critical fields and have been extracted reliably 
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[11, 12]. It is clear this time that we need a prolate order parameter to 
fit these data. We choose A(k) = - az 2 with a ~ 0.92 ~ 0.95. 

These values give a good fit of penetration depth data [11]. The magnetic 
penetration depths are related to the lower critical fields via the formulas 



H c cl {t) 



$0 



In 



Xg(t) 

at) 



In 



2ir\ a (t)\ c (t) V Ut)Ut) 



(1) 
(2) 



where X a (t) and X c (t) are the magnetic penetration depth with the super- 
current in the ab-plane and in parallel to the c-axis, respectively. They 
are related to the superfluid density via p s ,a{t) = A^(0) / A„(i) and p s , c (t) = 
X 2 c (0)/\l(t). Taking a = 0.95, which fits p a , a (t) from Ref. [11], we calcu- 
late p s ,c(t) an d obtain the temperature dependences of and from 
Eqs. (1) and (2), neglecting the temperature dependence of the logarithms. 
The result is shown in Fig. 1 along with the experimental results from Ref. 
[12]. From these fits we obtain H^(0) = 24 mT and H^(0) = 32 mT. 
From this the ratio of the relevant Fermi velocities can be estimated as 
Hci (0)/-^ci(0) ~ v c/ v a = 1-3. In other words, the corresponding Fermi 
surface is more isotropic and further v c > v a . This strongly suggests that 
the prolate order parameter we are considering has to be associated with 
the 7r-band. 




Fig. 1. Temperature dependence of the lower critical fields (dashed line) and 
(solid line) along with the corresponding experimental data from Ref. [12]. 
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4. Synthese 

We have seen so far that we need two energy gaps of different shape in 
order to describe H C 2(t,@) and H c i(t, 6). The temperature dependence of 
j(t) indicates that an oblate order parameter (A(/c) = l/Vl + az 2 ) dom- 
inates the behavior at high magnetic field. Also we need v c /v a = 0.48. 
This suggests the cylindrical Fermi surface associated with the cr-band as 
the carrier of this order parameter. On the other hand, for H c i(t,@) we 
need a prolate order parameter (A(k) = — az 2 ). Also, the anisotropy 

in H c i(t, 0) suggests the 7r-band as the carrier of this order parameter. In 
other words, if we assume that the high field properties are controlled by 
the oblate order parameter attached to the cr-band, while the low field prop- 
erties are due to the prolate order parameter attached to the 7r-band, we 
have a consistent picture for superconductivity in MgB2- 

We believe that this picture will be crucial to understand also a variety 
of anomalies observed in the vortex state in MgB2- 

We would like to thank M. Angst, F. Bouquet, S. Haas, A. Janossy, 
A. Junod, B. B. Jin and N. Klein for useful discussions on MgB2- 
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